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Time-dependent polarization relaxation behaviors induced by a depolarization field Ed were inves- 
tigated on high-quality ultrathin SrRuOa/BaTiOs/SrRuOa capacitors. The Ed values were deter- 
mined experimentally from an applied external field to stop the net polarization relaxation. These 
values agree with those from the electrostatic calculations, demonstrating that a large Ed inside the 
ultrathin ferroelectric layer could cause severe polarization relaxation. For numerous ferroelectric 
devices of capacitor configuration, this effect will set a stricter size limit than the critical thickness 
issue. 
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With recent breakthroughs in fabricating high-quality 
oxide films 0, 0- ultrathin ferroelectric (FE) films 
have attracted much attention from the scientific as well 
as application points of view. As the FE film thickness 
d approaches tens of unit cell length, the FE films of- 
ten show significantly different physical properties from 
those of bulk FE materials. Some extrinsic effects, es- 
pecially coming from FE film surfaces and/or interfaces 
with other materials, could be very important For 
some other cases, intrinsic physical quantities could play 
vital roles in determining the unique properties of ultra- 
thin films. 

Many FE-based electronic devices have the capacitor 
configuration, where a FE layer is inserted between two 
conducting electrodes. Then, polarization bound charges 
will be induced at the surfaces of the FE layer, but com- 
pensated by free charge carriers in the conducting elec- 
trodes. In real conducting electrodes, however, the com- 
pensating charges will be induced with a finite extent, 
called the screening length A. This will result in an in- 
complete compensation of the polarization charges. Such 
an incomplete charge compensation should induce a de- 
polarization field Ed inside the FE layer, with a direction 
opposite to that of the FE polarization P |{|. Therefore, 
Ed will appear in every FE capacitor, and its effects will 
becomes larger with the decrease of d [j| . (For a FE film 
without electrodes, there is no compensation for the po- 
larization bound charge, so the value of Ed will become 
even larger than that of the FE capacitor case.) Ed has 
been known to be important in determining the critical 
thickness and domain structure of ultrathin FE films 
, and reliability problems of numerous FE devices 

Recently, using a first principles calculation, Jun- 
quera and Ghosez investigated the critical thickness of 
BaTi0 3 (BTO) layers in SrRu0 3 (SRO)/BTO/SRO ca- 
pacitor Q. For calculations, they assumed that all of 
the BTO and SRO layers were fully strained with the 
SrTiO-3 substrate. By taking the real SRO/BTO inter- 
faces into account properly, they showed that Ed could 
make the ferroelectricity vanish for the BTO films thin- 
ner than 6 unit cells, i.e. 2.4 nm More recently, 



using pulsed laser deposition with a reflection high en- 
ergy electron diffraction monitoring system, we fabri- 
cated high-quality fully-strained SRO/BTO/SRO capac- 
itors on SrTi03 substrates with d between 5 and 30 nm 
With a very low leakage current, we could directly 
measure their P-E hysteresis loops Q. In this letter, we 
report the time-dependent polarization changes of the ul- 
trathin BTO films. We find that the net P of the ultra- 
thin BTO films decreases quite rapidly in time. We will 
show that the P relaxation should originate from Ed- By 
compensating for Ed with an external potential, we can 
determine the Ed values of the SRO/BTO/SRO capac- 
itors experimentally. These measured Ed values agree 
with the values from the electrostatic calculations. Fi- 
nally, we will discuss the effect of the P relaxation on a 
practical size limitation imposed on actual FE devices. 

In our earlier report Q, we obtained the thickness- 
dependent remnant polarization P r values from the P- 
E hysteresis loops, measured at 2 kHz in ultrathin FE 
films as thin as 5 ~ 30 nm. With further studies on the 
frequency dependence of the P r values in P-E hysteresis 
loops, as shown in Fig. 1(a) for a 15 nm thick BTO capac- 
itor, we found differences in the P r values when the mea- 
suring frequency is varied. These results suggest that the 
FE domain dynamics should play an important role for 
ultrathin FE films, where the FE domain wall motion is 
known to be strongly suppressed 13]. Note that the first 
principles calculation (FPC) and the Landau-Devonshire 
calculation (LDC) do not consider the domain dynamics, 
so their predicted polarization values should be called as 
spontaneous polarization P s . 

Since the P value significantly affects the subsequent 
analysis of P relaxation, precise determination of P s val- 
ues is necessary. To determine the precise values of P s , 
we applied pulse trains, which are schematically shown in 
the inset of Fig. 1(a) 14]. The interval between write and 
read pulses was set to 1 /zs to minimize the effects of the 
P relaxation, and the current responses under the read 
pulse were measured. The total amount of charge is ob- 
tained by integrating the current responses in time. The 
read pulses with different heights were used to obtain the 
linear part of the polarization under an external electric 
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field The P s values can be obtained by extrapolating the 
linear part of the polarization to zero electric field. The 
triangles (black) and circles (green) in Fig. 1 (b) show the 
P r values measured at 2 and fOO kHz, respectively. Also, 
the squares (red) show the P s values from the pulse test. 
The solid (green) and dashed (blue) curves show the the- 
oretical predictions from the FPC and the LDC 0, 
respectively, which take account of Ed- Note that neither 
of these theories can explain the thickness-dependence of 
P s quantitatively. However, it is known that the FPC 
predicts systematically somewhat lower bulk lattice con- 
stants compared to real values, so the compressive stress 
predicted by the FPC could be smaller than that in the 
fully strained sample, resulting in a smaller P s . To avoid 
this systematic error, we normalized the polarization val- 
ues to those of a 30 nm thick BTO capacitor. We found 
that the thickness-dependent scaling of P s also follows 
the FPC predictions quite well, as shown in Fig. 1(c). 

The large difference in P values between the 2 and the 
100 kHz tests indicates that there should be a strong 
change in the net P between 10 and 500 [is. Time- 



o Pulse Test 
/ / Linear fit - 




50 
40 
30 
20 
10 


( 

1.2 

_ 10 
E 

o 08 

CO 

& 06 

5T 

£ 04 
0.2 

0.0 






1 2 
E(MV/cm) 


„(b) 


■ 

■ ■ o - 


- / ■ 


A A 


. f 

1 B ° ^ 


^""^ LDC 

— FPC 
■ F s 

o P r (100 kHz)- 
A P^(2kHz) - 



10 20 
d (nm) 



30 



:(c) 






■ ■ 






- i m/ 






LDC - 




— FPC 




■ P 



10 20 

d (nm) 



30 



FIG. 1: (color-online) (a) Upper halves of hysteresis loops 
for 15 nm thick BTO capacitor at the measurement frequen- 
cies 2 and 100 kHz. The values of spontaneous polarization 
P a were obtained from linear extrapolation of pulse measure- 
ments, (b) The polarization values obtained from pulse mea- 
surements, first principles calculation (FPC), and Landau- 
Devonshire calculation (LDC). (c) Normalized behaviors of 
LDC, FPC, and P s to the values of BTO capacitors with 30 
nm thickness. 
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FIG. 2: (color-online) (a) Polarization relaxation in the 15 
nm thick BTO capacitor. The inset shows the schematics 
of measuring relaxation behaviors under an external voltage, 
(b) Slowing down of the relaxation behavior under external 
voltages in the 15 nm thickness BTO capacitor, (c) Thickness- 
dependent Ed in the ultrathin BTO capacitors. The Ed val- 
ues obtained from relaxation behaviors experimentally (red 
solid circles) and those from electrostatic calculations with the 
parameters determined and measured P s in this work (open 
squares). The (green) line indicates the Ed from electrostatic 
calculations with the polarization values obtained from first 
principles calculation. 



dependent P changes were investigated by applying two 
kinds of pulse trains, as shown in the inset of Fig. 2(a). 
For the write and the read pulses with the same (oppo- 
site) polarities, the amount of nonswitching (switching) 
P can be determined ^(J- The difference AP, between 
the switching and the nonswitching P should be twice 
as large as the net P. As shown in Fig. 2(a), AP de- 
creases quite rapidly for the film with d = 15 nm; AP 
falls to less than 10% of the P s value within a relaxation 
time t re i ax of 1000 s. As shown with the solid squares 
(black) in Fig. 2(b), AP decay follows a power-law de- 
pendence on t re i ax ■ Similar power-law decays of AP were 
observed for all the BTO films in the thickness range of 5 
~ 30 nm. Note that such a strong polarization relaxation 
could pose a serious problem in capacitor-type ultrathin 
FE devices. 

What is the origin of such strong polarization relax- 
ations? We thought that they could be closely related 
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to large Ed induced inside the BTO films. To verify 
this idea, we slowed down the relaxation phenomena by 
applying an external voltage, as shown in the inset of 
Fig. 2(a). The values of the applied external electric 
field E ext were obtained by dividing the applied external 
voltage by the corresponding film thickness. When the 
external field is applied in the opposite direction of Ed, 
the potential gradient inside the FE layer will decrease. 
Figure 2(b) shows that the slope of the power-law de- 
cay becomes smaller, as E ext increases. Assuming that 
the depth of the double-well potential for BTO ferro- 
electricity can be considered negligible compared to the 
effect of Ed, we approximately determined experimental 
Ed values from the applied electric field under which the 
slope becomes zero. Since Ed is proportional to P, the 
Ed value should increase slightly on application of E ext . 
After correcting this minor contribution, we could de- 
termine the Ed values, which are plotted as solid circles 
(red) in Fig. 2(c). 

From electrostatic calculations on the capacitor geom- 
etry, Mehta et al. showed that 



Ed = - 



P ( 2e F /d 
e e F \2e F /d 



e e /A 



(1) 



where d is the thickness of the FE layer, and e F and e e 
are the relative dielectric constants of the FE layer and 
the electrode, respectively To obtain theoretical E d 
values for our SRO/BTO/SRO capacitors, we have to 
know accurate values of e e , A, and e F . Unfortunately, 
the reported physical parameter values in the literature 
vary aaaaLJ Also, we could not find any definite 
experimental study on e e . 

To obtain the value of e e for an SRO electrode, we 
used optical spectroscopy. We measured the optical re- 
flectivity spectra of epitaxial SrRu03 films (thickness: 
about 0.5 /im) in a wide frequency region between 5 
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FIG. 3: (color-online) (a) Frequency-dependent dielectric 
functions of SRO. (b) Capacitance-electric field curve of the 
SRO/BTO/SRO capacitor with 5 nm thick BTO. 



meV and 30 eV and performed a Kramers-Kronig analy- 
sis to obtain the frequency-dependent dielectric function, 
e(uj) {=e'(uj) + ie"(ui)]. The details of these measure- 
ments and analysis were published elsewhere 0, H^| . 
The open squares in Fig. 3(a) and the inset show ex- 
perimental values of e'(w) and e"{uj), respectively. Note 
that e e in Eq. (1) represents the dielectric response from 
the bound charges, namely bound electrons and phonons. 
Since SRO is metallic, there should be a large contri- 
bution from the free Drude carriers, which masks the 
dielectric response from the bound charges. To obtain 
e e , we decompose e(w) into a free carrier contribution 
tcoherent(u) and a bound electron contribution tbound(u) 
by fitting the experimental e(u>) with a series of Lorentz 
oscillators, which are displayed as the dotted (blue) lines 
in the inset of Fig. 3(a). The dash-dotted (blue) lines in- 
dicate the bound electron contribution. From the dc limit 
of €bound(u), we could estimate that the bound electron 
contribution to e e is about 8.17. The phonon contribu- 
tion to e e was evaluated in a similar way by analyzing the 
phonon spectra and found to be about 0.28 Conse- 
quently, e e is determined to be about 8.45. 

Using the carrier density n = 1.2 x 10 22 /cm 3 of SRO 
|20| . the experimental value of e e , and the effective mass 
of an electron rripjt = 7m e , where m e is the mass of a 
free electron 0, [23] , we applied the free electron model 
and obtained A = 0.8 ± 0.1 A 0, H3- We also measured 
e F from the capacitance-electric field C-E curves of BTO 
capacitors. Figure 3(b) shows the C-E curve for the 5 
nm BTO capacitor. The C-E curve has the hysteretic 
behavior typical for a FE capacitor. The BTO capac- 
itors with 5 ~ 30 nm thickness show almost the same 
e F -E curves. The e F values can vary from 70 to 230 de- 
pending on the applied E. Since most of our experiments 
were performed under a finite applied field, which corre- 
sponds to a value between 1 and 2 MV/cm, the e F were 
estimated to be about 80 0. 

With the measured values of e e , A, and e F , we could es- 
timate the theoretical Ed values from Eq. (1) with the P s 
values obtained from the pulse test. The open squares in 
Fig. 2(c) are the theoretical Ed values. The solid (green) 
line shows the theoretical Ed values with the P s values, 
obtained from the FPC. These theoretical Ed values from 
the electrostatic model agree quite well with the exper- 
imental Ed values, determined from the polarization re- 
laxation. It should be noted that the Ed values are com- 
parable with or even larger than the measured coercive 
fields (in our samples, 300 ~ 400 kV/cm). These large Ed 
values can cause P reversal and FE domain formation, 
which will result in a reduction of the net P value as time 
elapses. The fact that two independent determinations 
provided nearly the same Ed values demonstrates that 
the polarization relaxation behavior should be dominated 
by Ed inside the FE layer. 

Note that the E^-induced AP decay comes intrinsi- 
cally from the incomplete compensation of the P charges 
(due to the finite screening length of the electrodes) in 
real conducting electrode, so that it will inevitably pose 
a fundamental limit for most FE device applications us- 
ing the capacitor configuration. This limitation should 
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be much more severe than that due to the critical thick- 
ness of the FE ultrathin films 0- Even if the FE film 
is thicker than the critical thickness, it is feasible that 
the -E^-induced AP decay is large enough to make the 
net P decrease significantly, resulting in retention failures 
for numerous FE devices. As d decreases, Ed increases 
significantly. With the current miniaturization trends in 
some FE devices, the large value of Ed should play a very 
important role in determining the ultimate size limits of 
FE devices. 

In order to reduce device failure due to the polariza- 
tion relaxation, we can try to select better electrode and 
FE materials. Noble metals, such as Pt, have been con- 
sidered better electrodes because they have high carrier 
density (resulting in A values smaller than that of SRO). 
However, the e e values of typical noble metals are much 
smaller than that of SRO, i.e. 8.45 |23, so Ed in capac- 
itors with noble metal electrodes can be large. For ex- 
ample, Ed in the range of 500 ~ 900 kV/cm is expected 
for a 15 nm thick BTO film with noble metal electrodes 
(typically, A = 0.4 ~ 0.5 A, e e = 2 ~ 4). Thus, the .Ed- 
induced P relaxation for the ultrathin BTO capacitors 
with the noble metal electrodes could be at least equal 
to or worse than that with SRO electrodes. Proper FE 
material selection can be another option. Since PbTi03 
is known to have a much deeper double-well potential 
than that of BTO [3111, 

the P relaxation should occur 
at a much lower rate even with the same value of Ed. 
Optimization of FE materials should be of great impor- 
tance for the improvement of ultrathin film nanoscale FE 
device performances. 

In summary, we demonstrated that the depolariza- 
tion field inside the ferroelectric film could cause a se- 
vere polarization relaxation. By slowing down the re- 
laxation under an external field, we could determine 
the depolarization field in a real capacitor of ultrathin 
SrRuOa/BaTiOa/Si'RuOa experimentally, which result is 
in good agreement with electrostatic calculations. Our 
investigation demonstrates that the depolarization field 
originates from intrinsic properties of electrode material 
such as the finite screening length and that the depo- 
larization field should play an important role in domain 
dynamics in ultrathin FE films. The polarization relax- 
ation due to the depolarization field could pose a serious 
size limitation for ultrathin ferroelectric devices. 
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